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We study zero-temperature hysteresis in random-field XY and Heisenberg models in the zero- 
frequency limit of a cyclic driving field. We consider three distributions of the random field and 
present exact solutions in the mean field limit. The results show a strong effect of the form of 
disorder on critical hysteresis as well as the shape of hysteresis loops. A discrepancy with an earlier 
study based on the renormalization group is resolved. 



I. INTRODUCTION 



Hysteresis is common in systems subjected to a cyclic force [l(. It means that the response to a changing force 

depends on the history of the force. In particular, the response in increasing force is different from that in decreasing 

force. This is caused by the delay in responding to the force. Theoretically hysteresis should disappear if the 

force changes sufficiently slowly but this often corresponds to unrcalistically long time periods. Several complex and 

disordered systems such as permanent magnets show hysteresis over the longest practical time scales. A microscopic 

Y^ • model of this phenomenon reveals critical points on the hysteresis loop where non-equilibrium susceptibility of the 

system diverges [2|, y, |5| . These points are characterized by a diverging correlation length, scaling of various quantities, 

(_J ■ and universality of critical exponents that is reminiscent of equilibrium critical phenomena. Universal behavior covers 

a wide class of materials, but aspects of universality are best examined in the framework of specific models. We focus 

on the zero temperature hysteresis in classical spin models in a quenched random field in the limit of zero frequency 

of the driving field. These restrictions are not so drastic as may appear at first sight. In a pioneering work, Sethna et 

al [2| used these simplifying features to examine magnetic hysteresis in the random field Ising model. Their model 

reproduces several experimentally observed features including the shape of hysteresis loops, Barkhausen noise, and 

return point memory [2|-|^] . They employ a Gaussian distribution of the random field with mean value zero and 

standard deviation a that plays the role of a tuning parameter in the model. They find a critical value a c such that 

_ ■, ' for a < er c , each half of the hysteresis loop has a first order jump in the magnetization at some applied field h. The 

size of the jump goes to zero as a — > a c from below. If h c is the critical field at which the jump vanishes, {h c , a c } is 

I/-) a non-equilibrium critical point showing scaling of thermodynamic functions and universality of critical exponents in 

^Q , its vicinity. There appears to be a fair amount of experimental support for the predictions of this model [3[ . Silvcira 

OA ' and Kardar [6[ generalized this model to n-componcnt continuous spins. Using the renormalization group approach 

they showed that the generalized model also shows critical hysteresis, and the exponents are independent of n above 

[~^. ■ 6 dimensions. The significance of 6 dimensions is that it is the upper critical dimension of the model with a Gaussian 

random field [7| . In 6 and higher dimensions the action is adequately described by a quadratic term, and higher order 

terms are irrelevant in the renormalization-group sense. The quadratic action can be solved exactly and predicts a 

second order transition as a — > a c from below. The quadratic action is n-dependent but the critical exponents are 

£> , not. 

Recently we solved [8| a slightly different variant of the n-component spin model exactly in the mean field ap- 
proximation based on infinitely weak but infinitely long range interactions. Our variant used randomly oriented but 
jrt ■ fixed length random fields and spins. We were surprised to find strikingly different results from those of Silvcira and 
Kardar [y] for their quadratic action. For n = 2 (XY model), we found a first order transition and no critical point. 
We also found wasp-waisted hysteresis loops in this case. For n — 3 (Heisenberg model) as well, we did not find a 
critical point of the familiar type but one with rather peculiar criticality. The question arises why the predictions of 
two qualitatively similar models are so different. Is it because of our variant of the random field or due to the use of 
soft (variable length) n-component spins in the field theoretic formulation of the renormalization group. We are now 
in a position to resolve this question. We recover the predictions of the quadratic action of Silvcira and Kardar if we 
keep the length of vector spins fixed but allow the magnitude of the random field to have a Gaussian distribution. 
We also consider a rectangular distribution of the magnitudes of random fields to examine if a distribution with a 
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compact support gives a different behavior than an unbounded distribution. Rectangular and Gaussian distributions 
are known to give qualitatively different results in the case of Ising spins [5, 9]. However, for XY and Heisenberg 
spins we find these two distributions yield similar results. 

II. EQUATIONS OF MOTION 

The model studied in [8| is characterized by the Hamiltonian, 

H = -jY,Si.Sj-Y,hi-Si-h.Y,Si (1) 

z,j i i 

Here {Si} are spins located at sites {i = 1,2,... N} of a d-dimensional lattice, {hi} are quenched random fields, and 
h is a uniform applied field; Si, hi, and h are n-component vectors with magnitudes \Si\ — 1, \h%\ = a%, and \h\ = h 
respectively. The vectors {hi} are randomly oriented but their magnitudes {<Zj} may have one of several different 
distributions. We consider three cases: (i) Oj is constant independent of the sites (a; = a), (ii) a, has a Gaussian 
distribution with standard deviation a centered at the origin, and (iii) a^ has a uniform distribution in the interval 
[0 < cii < A], and zero otherwise. The discrete time dynamics of fixed-length spins at zero temperature is given by, 

s t (t + i) = M\- / i (t) = J^4(i) + /T i + K (2) 

l/n*)l j 

Here fi{t) is the effective field at site-i at time t and the factor |/i(i)| in the denominator ensures that the dynamics 
does not alter the length of the spin vectors. We characterize the system by m(t) = jr^2ih.Si(t), the induced 
magnetization per spin along h. In the mean field limit a spin Si interacts with every other spin Sj with strength 
J = Jq/N. Let the applied field h be along the x-axis. Then the equations for the evolution of magnetization in XY 
and Heisenberg model in the three cases of the random field distribution mentioned above are [l0( : 

A. Randomly oriented fields of fixed magnitude cii — a. 

If the quenched fields {hi} are randomly oriented vectors of a fixed length a, the equations relating m(t + 1) to 
m(t) are, 

/. , ,, 1 f 2 * {J m{t) + h} + a cos a z 

m(t+l)= — / r dai XY model: n = 2 3) 

2ir J Q [a 2 + 2{ J m{t) + h}a cos a t + {J m(t) + h} 2 }^ 

, f\ 1 C { J()Tn(t) + h} + a cos ai . , rTT . n , . 

m[t + 1) = - / sin a,ida>i [Heisenberg model: n = 31 (4) 

2 Jo [a 2 + 2{J m(t) + ftjacosa, + {J m(t) + h} 2 ]? 

Equation (3) is a slight generalization of equation (9) in reference [8( for a = 1 . One may understand it intuitively 
as follows. The zero temperature iterative dynamics aligns Si(t + 1) along the effective field /»(£) at site-i. The 
effective field fi(t) is the sum of the mean field {Jom(t) + h} along the x-axis (the direction of the applied field) and 
a random field of magnitude a making an angle ai with the x-axis. Thus the components of fi(t) along the x and y 
axes are fo x = { Jom(t) + h} + a cos a^ and fi y = a sin ai. The component of the unit vector Si(t + 1) along the x-axis 
contributes to the magnetization m(t + 1). Its contribution is equal to the cosine of the angle that fi(t) makes with 



the cc-axis and therefore equal to fix/ \/ ff x + fi V - This explains the integrand in equation (3). The integral over ai 

amounts to taking an average over all sites to get the magnetization of the system. Equation (4) may be understood 
similarly. In this case the random fields as well as the spins are three component vectors and we need a polar and 
an azimuthal angle to specify their orientation. The polar angle a, of hi is measured from the x-axis which is again 
the direction of the applied field and the magnetization of the system. The azimuthal angle is integrated out trivially 
because the component of Si(t + 1) along the a;- axis does not depend on it leaving us with equation (4). 



B. Randomly oriented fields with a Gaussian distribution of a» 

Randomly oriented fields with a Gaussian distribution of at means that each Cartesian coordinate of the n- 
component field hi has a Gaussian distribution with average value zero and standard deviation a. The integrals 
in the equation relating m(t) to m(t + 1) are performed more conveniently in polar coordinates. We integrate over a^ 
with appropriate normalization and weight factors for two and three dimensional integrals respectively to get, 

n , -n 1 [°° -As a f 27T {Jom(t) + h} + a l cosa l 
m(t+l)=- — - e ^didai I — — — — ——rdai [XY: n = 2\ (5) 



27TCT 2 J Q J [ ? + 2{J m(t) + h} ai cosa; + {J m(t) + h} 2 }* 

f°° _<•? o, r {J m(t) + h} + a, cos a t 

/ e 2^ a j da,i I 

27rcr 3 Jo J a [a 2 + 2{ J m(t) + /i}a» cos on + { Jom(t) + h} 2 

C. Randomly oriented fields with a rectangular distribution of at 



{J m(t) + h} + a, cos a, , rTT . „. 

m(t+l) = / __ / e 2^ajaa; / — — ^— — — — — ■ - — ; suiajdaj [Hcisenberg: n = 3J 

(6) 



The equations for this case are obtained on the same lines as in the preceding section by integrating over a, with a 
uniform distribution in the interval [0 < ai < A] . 



1 f A , f 27r {J m(t) + h} + a t cos a t 

— pr / CLiddi I 

7i" A Jo Jo [af + 2{J m(t) + h}a t cos a t + { J Q m(t) + h} 2 ] 



m(t + 1) = rT5 / **»* / f ,.. fT ^^^JX^^_ __ ^ [XY: n = 2] (7) 



/ \ 3 f n , f* {Jom(t) + h} + o, cos a, , , TT . . „. .„. 

m(t + l)= —it / ajdai / — sin ondan [Hcisenberg: n = 3J (8) 



2A 3 7 ! 7 [ a 2 + 2{ J m(t) + ft,}a, cos a, + {J m{t) + h} 2 

III. HYSTERESIS IN XY MODEL 

Although the integral in equation (3) can be written in terms of elliptic functions of the first and second kind, it has 
an appealing geometrical interpretation @,[ll( in its present form. The left- hand- side represents the average projection 
of a spin along the cc-axis. This is proportional to two forces acting along the x-axis (i) mean field { Jam(t) + h} and 
(ii) a random field "acosai". These account for the numerator in the integrand. The denominator ensures that the 
resultant of the mean field and the random field vectors is a unit vector. The orientation of the random field changes 
from site to site and the integral over on represents an average over it. In order to examine the structure of a fixed 
point to* of equation (3), it is convenient to set Jo — 1 an d ft = so that to* is equal to the mean field. Consider 
two circles of radii unity and a respectively (a < 1) with their centers separated by m* on the £-axis as shown in 
figure (1). The spin lies on the larger circle and the random field on the smaller one. The spin making an angle 9 
with the z-axis generally cuts the smaller circle at two points. Correspondingly two orientations of the random field, 
one making an angle a and the other a + (3 with the a;-axis produce the same magnetization to = cos 9. Using the 
properties of triangles in figure (1), we get j3 + 2(q — 9) — -k and sin(7r — a + 9)/m* = s'm9/a. These two relations 
can be combined to give, 

to 2 = 1-f— Ycos 2 ^- (9) 

\m*J 2 v ; 

In the limiting case j3 — when the spin vector is tangential to the random field circle, to approaches its minimum 
value m m i n . For a given a, the fixed point to* is an average over various values of m ranging from to = m m i n (for 
f3 = 0) to m = 1 corresponding to a = and /3 = ir. Note that the integrand in equation (3) has the same value for 
on = 7r as for on = 0. Figure (2) shows the fixed points of equation (3) with Jo = 1 and h = for increasing and 
decreasing a starting from to* = 1 at a = and changing a in small steps 8a. At each step, the fixed point at the 
preceding step is used as a starting point and a is held fixed during the iterative dynamics till a new fixed point is 
reached. We see that to* decreases with increasing a; at a w 0.67, to* jumps down from to* w 0.74 to to* = and 




FIG. 1: (Color online) The figure shows a geometrical relationship between the spin, the random field, and the mean field at 
the threshold of instability in increasing a for Jo = 1 and h — (see text). 



stays zero thereafter. On decreasing a, m* stays zero for a < 0.5, jumps up to m* sa 0.92 at a = 0.5 and increases 
along the old trajectory as a — > 0. The region 0.5 < a < 0.671 is bistable. It contains two lines of stable fixed points 
separated by a line of unstable fixed points. The unstable fixed points separate the domains of the stable fixed points. 
In the bistable region, the system is in the domain of m* ^ in increasing a, and m* = in decreasing a. This is the 
basic mechanism of hysteresis in our model. The boundaries of the bistable region can be understood geometrically. 
As a increases from zero the random field circle increases and its center shifts closer to the center of the larger circle. 
The angle a corresponding to m* also increases until a = tt/2 at m* ~ 0.74 at a « 0.67. This is depicted in figure (1). 
For larger a, m* = and the two circles are concentric in the stable state. The first order jump in the magnetization 
at {h = 0, a w 0.67} gives rise to hysteresis loops in cyclic fields as shown in figure (3). Note that m* — at h = 
if a > 0.67 but this does not mean the absence of hysteresis but rather wasp-waisted hysteresis loops as explained in 
reference Q . 

If the random field at each site has a random orientation at and a Gaussian distributed magnitude a% with mean 
zero and standard deviation a, our numerical work suggests that the wasp-waisted hysteresis loops disappear and the 
first order jumps in the magnetization gradually diminish with increasing a and vanish as a — > a c from below. We 
find cr c « 0.62. Figure (4) shows the hysteresis loops for two values of <r, (i) a < a c and (ii) a > a c . The critical 
behavior may be obtained analytically. It is again useful to set Jo = 1 and h = in equation (5), and then expand it 
in powers of m(t) in the limit m(t) — > 0. We split the range of the integral over a» in two parts (i) zero to m(t) and 
(ii) m(t) to oo and extract the leading terms from both intervals. Thus, 
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(10) 



After performing the Gaussian integrals and simplifying, we get 



m(t+ 1) 



m(i) 



i(t) 



(11) 



This recursion relation shows that a has a critical value a c 

\l/2 



0.6267; 



= if 



<y > (Jr. 



and 



v / 8cr c ((T c — a) * as a — > a c from below. Thus there is a second order phase transition at a — a c . These predictions 
are born out by numerical solution of the equation as shown in figure (5). 

For randomly oriented fixed- length fields {a,; = a}, the corresponding equation for the magnetization is, 



m(t+ 1) 



m(t) 



1 (m(t) 
8 1 a 



(12) 



As explained in reference [8(, the positive sign of the cubic term in the above equation rules out a continuous 
transition from in* = 0. There is a first order transition in this case. One may ask if the difference in hysteresis 
for a Delta function distribution {a^ = a} and a Gaussian distribution comes from the difference between a sharply 
localized distribution and an unbounded distribution. In order to examine this question we first relax the equality 
a, = a and allow a, to be uniformly distributed in a narrow range [a — e to a + e] . We find that the hysteresis for 
a narrow rectangular distribution around a, = a is qualitatively the same as that for a, = a in the limit e — > 0. In 
the same vein, we also consider the distribution at = ■% if < a^ < A, and ai = otherwise. For this distribution, 
equation (7) takes the following form in the limit m(t) — > 0, 



m(t+l) 



m(t) 



1 
1- - 



m{t) 



(13) 



The above recursion relation has a critical point at A c = 1; m* — if A > A c , and m* = y / 8A 2 (A c — A) as 
A approaches A c from below. Figure (6) shows m*(h = 0) as A increases from zero to A > A c . Thus the critical 
behavior of the model with a uniform bounded distribution of {ai} is the same as for a Gaussian distribution. 



IV. HYSTERESIS IN HEISENBERG MODEL 



We now examine hysteresis in the Hciscnbcrg model with each of the three distributions considered above, 
random field vectors of a fixed length \dl\ = a and A(i) = Jom(t) + h, equation (4) simplifies to the following: 



For 



0.75 



0.5 - 



0.25 




FIG. 2: (Color online) Magnetization curves in increasing and decreasing a at h = 0. The thick dark (blue) line shows the 
magnetization in increasing a. Starting with m = 1 at a — 0, it decreases with increasing a, drops to zero at a » 0.67, and 
stays zero for larger values of a. For decreasing a, it remains zero for a < 0.5 along the thick gray (green) line, jumps up on the 
thick dark (blue) curve at a = 0.5 and follows it till the end. The thin dark curve (red) shows a line of unstable fixed points 
in the bistable region 0.5 < a < 0.671. On decreasing a, the starting magnetization is always below the line of unstable fixed 
points, and therefore in the domain of the m* = fixed point. 
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FIG. 3: (Color online) Hysteresis loops for the XY model with randomly oriented fields of length a=0.5 (dark/red line), 0.65 
(thicker gray/green line), and 0.9 (thicker dark/blue line). At a = a c (~ 0.67), \m(h = 0)| jumps from m ~ 0.74 tom = (see 
figure 2). Therefore hysteresis for a = 0.9 has a wasp-waisted loop with zero hysteresis in small applied fields around h — 0. 



m(t+l) = ^^ if \A(t)\ <a 
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3\A(t) 



if A(t) > a 



if A(t) < -a 



(14) 



Setting J = 1 and h = 0, the fixed point in the small m limit is determined by the recursion relation m(t + 1) = 
2m(£)/3a. Thus m* = is stable fixed point if a > 2/3. At a = 2/3 any value in the range — I < m* < | satisfies the 
fixed point equation at h = 0. This is a peculiarity of mean field hysteresis in the Heisenberg model with randomly 
oriented fields of fixed length a. If a > |, |m*| > | and approaches unity as a — >• 0. Therefore as /i is cycled between 
— oo and oo, we get hysteresis if a < 2/3, no hysteresis if a > 2/3, but hysteresis docs not vanish at a = 2/3 in any 
familiar fashion of first or second order transition. 

If the magnitude of the random field has a Gaussian distribution, the appropriate recursion relation in the small m 
limit is obtained from equation (6). Setting Jq = 1 and h = 0, we get 



m(t+l) = J~ 

V IT <J a 



<l. 



m _ „f r l 

e ^ajddi { 1 - - . 

3 \mi(t) 




' afdcLi 



2m % {t) 
3ai 



(15) 



This simplifies to, 



m(t + l) 



m(t) 



a 



1 (m(t) 
10 1 a 



+ 



(16) 



In contrast to equation (14), the above equation shows a square root singularity in critical hysteresis at a c = 
y/8/9n sa 0.5319; m* = if a > a c ; and m* = \/10i7 c (a c — a) 1 / 2 as a — > a c from below. Therefore magnetization 
curves for a > a c should be reversible and those for a < a c should show hysteresis loops. Figure (4) shows the results 
for a = 0.5 and a = 0.7 respectively. Figure (5) shows the critical behavior of m*(h = 0) as a — > a c . The fit with the 
expression derived above is good over a rather wide region a < a c . For comparison, we have also shown the result 
for the Gaussian random field Ising model in figure (5). In the mean field theory, the magnetization vanishes with 
a square root singularity as a — > a c irrespective of n (n = 1, 2, 3) but a c decreases with increasing n. This may be 
expected because spins with larger number of components have more freedom to disorder. The Gaussian random field 
Ising model has been investigated extensively analytically as well as numerically. Its predictions are in reasonable 
agreement with the experimental observations. It is argued [6| that even for XY and Heisenberg models, if the 
critical hysteresis occurs at m* ^ or h ^ 0, it may show critical exponents of the Ising model because the non-zero 
magnetization or applied field picks a unique direction in the system. It is not clear why the Gaussian distribution 
should be best suited for comparison with experiments, but there is comparatively little study of other distributions 
of the random field. A numerical study shows that critical exponents of three dimensional random field Ising model 
with a Gaussian distribution are significantly different from those of the same model with a bimodal distribution of 
random fields |12| . 

In the case of a rectangular distribution of a*, the equations are easily integrated. We get, 



m(t + 1) = -rm(t) 



A 



1 
1- - 

5 



m(t) 
~A~ 

A 
m(t) 



if m(t) < A 



if m(t) > A 



(17) 



These equations show that similar to the case of XY model with a rectangular d istribution of a,, the Heisenberg 
model shows a critical point at A c = 1. The fixed point is m* = y/5(A c — A)A 2 if m* < A giving a square root 
singularity similar to that for a Gaussian distribution. If m* > A, m* is determined by the real stable root of the 



cubic equation to 3 
zero to A > A c . 
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0. Figure (6) shows the fixed point magnetization at h = as A is increased from 
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FIG. 4: (Color online) Hysteresis in XY and Heisenberg models with randomly oriented fields of Gaussian magnitude. Results 
are shown for two values of a in each case; a = 0.5 (a < a c ) and a = 0.7 (a > <t c ). The gray (green) curves show the hysteresis 
loop and the reversible magnetization for the XY model for a = 0.5 and a = 0.7 respectively. The dark (blue) curves show 
similar results for the Heisenberg model. The Heisenberg model has a smaller hysteresis loop for a — 0.5 and its reversible 
magnetization for a — 0.7 has a smaller slope at h = 0. 
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FIG. 5: (Color online) Fixed point magnetization m*(h = 0) vs a starting with m = 1: Ising model (the curve to the 
right/blue; a c = a/7t/2 ~ 0.8), XY model (the middle curve/green; a c — \/n/8 « 0.6267), and Heisenberg model (the curve to 
the left/red: a c = ^ys/9^v ~ 0.5319). The expression (a c — a) 1 ' 2 with prefactors \Z6a c , \/8a c , y'lQtJc (see text) is superimposed 
on each curve respectively. The square root singularity is seen to fit the data over a rather wide critical region. 
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FIG. 6: (Color online) Fixed point magnetization m*(h = 0) vs A starting with m = 1: XY model (the gray/green curve; 
A c = 1), and Heisenberg model (dark/blue curve: A c = 1). The expression \/8A(A c — A) 1 ' 2 is superimposed on the XY result 



to bring out the square root singularity in the critical region. For the Heisenberg model m* = Ay5A 2 (A c 
exactly for m* < A (see text). 
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V. CONCLUDING REMARKS 

This work has been motivated by the need to resolve the discrepancy between two exactly solved models of critical 
hysteresis in n-component spins (n = 2, 3) placed in a quenched random field. A calculation [6| based on the RG 
approach predicts a second order transition with n-independent critical exponents if d > 6. This calculation uses a 
Gaussian distribution of the random field. The second calculation Q does not involve d explicitly but considers the 
mean field approach of each spin interacting equally with every other spin. It employs randomly oriented unit vector 
fields and predicts a first order transition and wasp-waisted loops if n = 2, and a peculiar transition if n = 3. We have 
shown that if the random fields in the second variant of the model have a Gaussian distribution, the mean field theory 
reproduces the results of the RG approach. This is reassuring and may have been anticipated. However an explicit 
verification is valuable because the two models and methods of solution are not identical. The RG approach uses 
variable length spins with nearest neighbor interactions on a d-dimensional lattice while the spins in the mean field 
approach are unit vectors. The striking difference in hysteresis between Delta function and Gaussian distribution of 
the magnitude of random fields is not generic to a localized and an unbounded distribution. We find that a rectangular 
distribution with a compact support yields the same critical behavior as the Gaussian distribution in the mean field 
theory. Thus the calculation presented here restores some expectations of the universality of non-equilibrium critical 
behavior but it also raises other questions that we discuss in the following. It brings out an unexpectedly strong 
dependence of the shape of hysteresis loops and critical hysteresis on the distribution of the quenched field in the 
system at the mean field level. These results have bearings on hysteresis in magnetic systems beyond the mean field 
theory [2, 0, mM, fl3H17l | as well as other adiabatically driven disordered systems |18l - [26l |. We hope the analysis 
presented here may contribute to a better understanding of hysteresis experiments and also motivate future work to 
explore what features of quenched disorder are crucial in determining the observed phenomena. 

It would be interesting to understand our analytic results on general grounds such as the symmetry of the system and 
universality classes of critical behavior. Why is it that randomly oriented fixed-length fields produce different critical 
behavior than randomly oriented fields with Gaussian distributed magnitudes? Why is it that Gaussian distributed 
magnitudes produce the same critical behavior as uniformly distributed magnitudes? We know that critical hysteresis 
for Ising spins is different for a uniform bounded distribution of fields than for a Gaussian distribution [5|, [9| . An 
exact solution |5[ of the model on a Bethe lattice with coordination number z reveals that there are no jumps in the 
magnetization nor there is a critical point on the hysteresis loop if the fields have a Gaussian distribution and z < 3. 
For a uniform bounded distribution of fields there arc jumps in the magnetization even for z = 3, and the jumps go to 
zero discontinuously as the width of the distribution is increased. This is due to instabilities in the equations of motion 
at large quenched fields on the boundary of the distribution. There is also a renormalization group argument [J, |9( 
that for Ising spins the crucial feature of the probability distribution of random field is its second derivative at zero 
field. In this formalism the difference in the critical behavior of the models with uniform and Gaussian distributions 
is attributed to the difference in the second derivative of the probability distribution of the random field at zero 
field. Although the renormalization group analysis as well as an explicit solution of equations of motion with a given 
initial condition lead to similar conclusions for critical hysteresis in the case of Ising spins but it does not suggest 
any simple renormalization group argument why uniform and Gaussian distributions should yield the same critical 
hysteresis in XY and Heisenberg models. The results presented here are based on an exact solution of the equations of 
motion in the mean field approximation. The renormalization group tackles the problem rather indirectly. It recasts 
the equations of motion into a path integral that is a sum over all paths of the exponential of an effective action. 
The effective action carries extra baggage compared with the equations of motion. This baggage (extra terms and 
parameters) comes from the use of soft spins and auxiliary fields needed in the exponentiation of delta functions. It 
is expected that the extra terms are irrelevant in the renormalization group sense and do not influence the critical 
behavior of the model. However it is not clear to us at this time if the difference between a uniform and Gaussian 
distribution of quenched fields is irrelevant for spins and fields having more than one component. 

An intuitive perspective of results presented in this paper may be obtained by paraphrasing them in terms of 
energy barriers in the dynamics of n-component spins in quenched fields comprising randomly oriented ?i-component 
(n > 2) vectors of length a. Specifically, consider the dynamics of the system in the absence of an applied field for 
two different initial conditions (i) all spins pointing along the negative x-axis (mo = — 1), and (ii) all spins pointing 
along the positive rr-axis (mo = 1). In each case the relaxation dynamics takes the system to a stable fixed point 
whose magnetization is close to the initial magnetization if a << J where J is the ferromagnetic exchange coupling 
between spins. The two fixed points are separated by a large gap in their magnetizations and consequently by a 
large energy barrier between them. As a is gradually increased (keeping J fixed) the gap between the two fixed 
points decreases. This is understandable because larger a means larger disorder in the system and the fixed point is 
expected to move farther away from the initial uniform state. In the mean field theory of XY spins, the gap goes 
to zero discontinuously at a critical value of a = a c . If the lengths of the random field vectors {a^} have a uniform 
distribution in the range < a, < A or a Gaussian distribution in the range < a, < oo with standard deviation u, 
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the gap goes to zero continuously at a critical value of A or a as the width of the distribution is gradually increased. 
This too is understandable because the distributed magnitudes, particularly very small magnitudes in the vicinity of 
zero magnitude would lower the energy barriers for rotation of spins at their sites. The barrier for rotation of a spin 
depends on the magnitude and direction of the quenched field at its site as well as the orientation of its neighbors. 
A connected cluster of low at sites is likely to have a much lower barrier for collective movement in an avalanche 
than if all a% were equal to each other and of the order of J. We may also expect the difference between uniform 
and Gaussian distributions of {a^} to be less significant for XY spins than for Ising spins that can only flip but not 
rotate. Similar considerations apply for Hciscnbcrg spins as well where the distributed magnitudes {a^} would push 
the energy barriers even lower due to an additional degree of freedom for rotation. In a future study we plan to 
investigate the extent to which the predictions of the mean field theory may apply to hysteresis in random field XY 
and Hcisenberg models with short range interactions on periodic lattices. 
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